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1 Abstract-Nike Apache 14.439 was launched from Wallops Island at 0003 EST on 

| 1 November 1972, a very disturbed night (K =8). A Geiger counter in 
I P 

the payload detected electrons (>70 keV) with a maximum flux of 1086 + 

i -2 -1 -1 

I 26 cm sec ster . The height-averaged ionization rate in the upper E region 

is calculated from the measured electron density profile and has a value of 
-3 -1 _ . 

35 cm sec . This ionization rate can be' reconciled with the observed flux 

of electrons O?0 keV) if the spectrum (>2 keV) is of the form J(>E) - J 

■ 

j exp(-E/E Q ) with E q equal to 8.3 keV. The ionization rate on this and other 
I nights is found to be strongly dependent on geomagnetic activity. It is 
suggested that energetic electrons are the principal source of ionization at 
midlatitudes in the upper E region near midnight, even under rather quiet 
geomagnetic conditions. 

1. Introduction 

Rocket observations at Wallops Island (38°N, 75°W; geomagnetic latitude 49°N) 
have shown the presence of a layer in the upper E region near midnight 
(Cartwright, 1964; Smith, 1970). This intermediate layer is formed by the 
vertical redistribution of ionization under the influence of the neutral winds 



(Const antinides and Bedinger, 1971; Fujitaka et. al., 1971; Chen and Harris, 
1971). The source of ionization has been presumed to be solar UV radiation 
resonately scattered by the geocorona, principally Lyman-a (1216 A), which 
ionizes nitric oxide (Swider, 1965) and Lyman-0 (1026 A), which ionizes 
molecular oxygen (Ogawa and Tohmatsu, 1966) . 

The possibility exists, particularly on disturbed nights in the sub- 
auroral zone (geomagnetic latitude 45° to 60°), that energetic electrons 
provide an additional source of ionization in the upper E region and can 
therefore influence the intermediate layer (Manson, 1971). The rocket experi- 
ment which is described here was undertaken to examine the role of energetic 
electrons as a source of ionization in the midlatitude nighttime ionosphere. 

The rocket launch time was chosen to be midnight to minimize the effect of 

UV sources of ionization and the launch held for a night of considerable 
geomagnetic activity. 

After presenting the observations from this rocket launch, we consider 
the interpretation of the observed enhancement in electron density and deduce 
a flux of electrons with energies in the range 2 to 10 keV necessary to explain 

the effect. Tins flux is shown to be reasonable for the disturbed conditions 
of that night. 

In a later section of the paper we consider previous observations of the 
intermediate layer near midnight. Making proper allowance for the vertical 
redistribution of ionization resulting from the effect of the horizontal 
neutral wind, we deduce the average ionization rate of the upper E region. 

Tins is found to be well correllated with the planetary index of geomagnetic 
activity, which leads us to suggest that, near local midnight, energetic 
electrons are an important source of ionization in the upper E region, not 
onj)' on disturbed nights but also for quiet conditions. 
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2. Instrumentation 


Nilce Apache 14.439, instrumented to obtain the electron density profile and 
to measure the flux of energetic electrons, was. launched from Wallops Island, 
Virginia, at 0003 EST (0503 UT) on 1 November 1972, during a magnetic storm. 
The value of for the three-hour interval (03-06 UT) including the launch 
time was 8. The rocket achieved an apogee of 195.6 km and, in mid-flight 
had a spin rate of 7.3 rps. 

The payload included a Langmuir probe using the nose tip electrode 
(Smith, 1969). This is operated at a fixed potential of +4.05 volt until 
100 sec after launch (127 km altitude) . For the remainder of the flight the 
probe alternates between the fixed-voltage mode of 1.5 sec duration and a 
swept mode of 0.5 sec duration. The constant of proportionality between probe 
current in the fixed- voltage mode and electron density changes with altitude 
and varies from flight to flight. The probe is therefore calibrated in flight 
by a propagation experiment. For this flight two frequencies were used: 

2225 kHz and 3385 kHz. The propagation experiment at each frequency uses 
differential absorption and Faraday rotation data (Mechtly and Smith, 1970) . 

The payload also included a Geiger counter with a window of beryllium, 

2 3 

area 0.203 cm and thickness 5.1 x 10 cm. The counter was originally designed 
for measurement of the daytime flux of solar X-rays in the wave length range 
1 to 8 A (Accardo et. al., 1972). It is also sensitive to electrons with 
energies greater than 70 keV, for which the geometric factor is 0.303 cm 2 ster. 
Although such relatively high-energy electrons are important only in the D 
region the counter was included in the payload as a direct indicator of a 
particle precipitation event. The window of the Geiger counter is oriented 
perpendicular to the spin axis of the rocket and is covered, during the launch 
phase, by a door carrying a weak radioactive source. The door and source are 
ejected 41 sec after launch (43 km altitude) . 
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3. Observations 


3.1 Flux of energetic electrons 

The count rate measured during the flight is shown in Figure 1. Each 
data point is the average for a five-second interval. The uncertainty in 
count rate is indicated by error bars at representative points. The data 
below 60 km on the descending trajectory are of poor quality and are not 
included in the figure. 

The count rate of the detector at sea-level (before launch) is 3.3 sec -1 , 
of which the source on the door contributes 3.0 sec" and the cosmic-ray 
background 0.3 sec . The dashed line in Figure 1 connects points which have 
been corrected for the radioactive source carried in the early part of the 
flight. The resulting variation with altitude up to about 80 km, including 
the Pfotzer-Regener maximum near 20 km, is consistant with previous observa- 
tions of cosmic ray fluxes at midlatitudes (Van Allen and Tatel, 1948). In 
other (daytime) flights with identical counters the background count rate, 
due to cosmic rays, has been found to be about 10 sec at altitudes greater 
than 40 km. 

The excess counts on this flight are attributed to electrons with energies 

greater than 70 keV . The flux increases with increasing altitude in the range 

80 to 140 km both on ascent and descent. While the rocket is above 140 km 

tile re is no obvious relation with altitude, the flux being variable and 

gradually decreasing. The count rate reaches a maximum value of 339 ± 8 sec” 1 , 

which, when corrected for cosmic rays, gives a flux (>70 keV) , assumed isotropic, 

of 1086 ± 26 cm “ sec ster" . In the period that the rocket is above 140 km 

“ ^ "1 1 

tht; minimum flux is 132 + 9 cm sec ster . These fluxes have been compared 
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with other midlatitude data collected by Potemra and Zmuda (1970) and are 

found to be close to their upper limit, confirming that this occasion was 
very disturbed. 

3.2 Electron Density 

The electron density profile from this flight is shown in Figure 2, 

together with that from Nike Apache 14.394, also launched at Wallops Island 

near midnight, but on a much quieter night (K : 2+). The earlier profile is 

J: 

typical of previous observations of the E region near midnight. The main 
features of these are the lower irregular region between 90 and 110 km, the 
valleys at about 120 and 180 km, with the intermediate layer between. The 
profile from the disturbed night shows a generally greater electron density 
except between 93 and 105 km with much irregular structure present in the 
upper part of the profile. There is some indication of the intermediate 
layer but its presence is not as obytous a.s in any previous profile. 

It is clear that the enhancement below 93 km is caused by energetic 
electrons. In the subsequent sections of this paper we shall examine the 
possibility that energetic electrons also are responsible for the enhancement 
above 105 km. That there is no substantial enhancement between these two 
altitudes is tentatively explained by the presence of a source of ionization 
from meteors, which predominates in the range 93 to 105 km. The present paper 
is, however, limited to detailed consideration of the region above 105 km. 

4. Analysis 

4.1 Introduction 

Neutral winds plan an important role in determining the electron density 
profile of the nighttime E region. In this section, we illustrate the forma- 
tion of the intermediate layer by solving the continuity equation for an 
assumed wind structure. Then we show that, by suitable manipulation of the 



continuity equation, we may perforin the inverse procedure and obtain the 
ionization rate, without assuming a value for the vertical ion drift velocity. 

4.2 Effects of winds in the E region 

Tine steady-state continuity equation for electrons in the E region is: 


q - 



_d 

dz 


(Nw) 


dz 



0 


Cl) 


in which q is the ionization rate, a is the recombination coefficient, N is 
the electron density, w is the vertical ion drift velocity, D is the ambipalar 
diffusion coefficient and z is the altitude. The small gravitational 
diffusion term has been neglected. 

The vertical ion drift velocity is produced by the neutral wind and can 
be calculated using the formulation of MacLeod (1966). Electric fields may 
also contribute to the vertical ion drift velocity; the subsequent analysis 
is independent of the mechanism producing vertical ion drift, however. 

For purposes of illustrating the effect of winds in the upper E region 
we have taken a model in which the vertical ion drift velocity, at altitudes 
between 120 and 180 km, is given by 

w = 20sin[7r(150 - z)/30] (2) 

-1 

where w is in m sec and z is in km. This forces the peak to be located at 
an altitude of 150 km and the valleys to be located at altitudes of 120 kin 
and ISO km. The choice of a maximum vertical velocity of 20 m sec", is based 
on the work of Gonstantinidcs and Bedinger (1971). 
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The steady-state continuity equation (1) has been solved numerically 

using the Runge-Kutta technique. For simplicity both a and D have been given 

-7 3 - 1 

values appropriate to 150 km: a = 1,90 x 10 cm sec (Fujitaka et al . , 

8 9 _ 2 

1971) and D = 2 x 10 cm sec . The latter is twice the molecular diffusion 
coefficient given by Golomb and MacLeod (1966) . 

The solution is shown in Figure 3 for five values of ionization rate: 

(1) 0.4, (2) 4.0, (3) 40, (4) 400 and (5) 4000 cnf 3 sec -1 . The effect of 
omitting the diffusion term in the continuity equation is represented by the 
thin line. The full solution, with diffusion, is shown by the thick line. 

The difference is most noticable for the smaller production rates, as on 
quiet nights, and is completely negligible for the largest production rate, 
which is representative of daytime conditions. It is important to note that 
diffusion cannot be neglected in analysis of the intermediate layer at night. 

This figure confirms that the wind becomes less effective in determining 
the electron density profile as the production rate increases, explaining, at 
least qualitatively, why the intermediate layer was less obvious on the 
disturbed night. 

Profiles of electron density, such as those shown in Figure 3, may be 
analyzed to obtain the ionization rate. This has previously involved a tedious 
modelling procedure. The new method described here follows from integration 
of the continuity equation between suitably chosen limits of altitude. It does 
not require independent knowledge of the vertical ion drift velocity. The 
ionization rate is determined as the height -integrated value appropriate to the 
range of altitude. 
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4.3 Calculation of ionization rate 

Consider the steady-state continuity equation (1) integrated between the 


altitude limits and z ? : 


f 22 

(q-aN^)dz - 

* 1 
^ / 

Z 2 ' 

+ 

z i 

L j 

1 Z ]_ 


1 Z 2 
,dN ! 1 


tlz 


= 0 


(3) 


Reference to the calculated profiles of Figure 3 shows that the extrema 
of electron density (i.e. ^ = 0) occur at nodes of the vertical ion drift 
velocity (i.e. w = 0). Thus, if the limits of integration are altitudes 
z^ and z^, which are extrema of electron density and nodes in the vertical 
ion drift velocity, the terms in square brackets in equation (3) are zero at 
both limits. Therefore: 


q dz = 


aN 2 dz 


(4) 


1 1 

Thus, the height-integrated ionization rate in the layer is obtained 
without assuming any model for the vertical ion drift velocity. The proper 
choice of limits of integration has enabled us to average out the effect of 
vertical ion drift (whether caused by neutral winds or electric field, or 
both) and diffusion; they have not been neglected. 

In the absence of further information or assumption about the variation 
of ionization rate with altitude it is convenient to define an average ioniza- 
tion rate given by: 


z 

( 2 


q = 


(Z 2 - Z P j 


(cdT) dz 


(5) 


S 



Ihe conditions for the limits of integration are met at both lower and 
upper valleys and at the peak of the layer. Thus the average ionization rate 
may be computed for (1) the layer below the peak, (2) the layer above the 
peak, and (3) the whole layer. 

Further discussion of this method of analysis and its extension to obtain 
the vertical ion drift velocity are given in Geller et_ al. (1974). 

5. Discussion 

5.1 Ionization rate on the disturbed night 

At night, in the altitude range of interest here, the pi'incipal ion is 

NO (Holmes et^a^., 1965). The loss process is dissociative recombination of 

N0 + , the coefficient (a) being given by 4.5 x 1(T 7 (300/T ) cm 3 sec” 1 

(Biondi, 1969). The electron temperature (T ) above 120 km varies with altitude 

and geomagnetic activity. For purposes of calculation we have used data from 

Evans (1973) at Millstone Hill and find that the nighttime (2100-0300 EST) 

electron temperature at 200 km can be well represented by (112 K + 588) °K 

for eight occasions in 1966 and 1967 with a range, of values of of 1- to 5-. 

The correlation coefficient is 0.90. Following Evans we assume T - 355°K at 

e 

120 km. Then a linear interpolation between these two altitudes gives 

T (z) = (z - 120) x (1.40 K + 2.91) + 355°K. (6) 

e p 

Using Biondi f s formula for a and allowing the extrapolation of electron 
temperature to K = 8 we obtain a as a function of altitude and geomagnetic 

Jr 

index. 

This value of a has been used in integrating the electron density profile 

according to equation (5) to obtain the ionization rate. For the profile of 

1 November 1972 between the altitudes 115 and 175 km, the ionization rate is 
-3 -1 

found to be 35 cm sec 
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5.2 Electron energy spectrum 

Calculations such as those by Rees (1963) show that electrons with energies 

in the range of 2 to 10 keV are the most efficient as a source of ionization in 

the upper E region, at about 140 km. In this range of energy the ionization 

-5 -3 -1 

rate for unit flux is about 1 x 10 cm sec . Thus to produce an 
... -3-1 

ionization rate of 35 cm sec requires a flux of electrons, with energies 

6 2 X "1 

between 2 and 10 keV, of about 3.5 x 10 cm* sec" ster" . 

The energy spectrum of electrons between 2 and 70 keV was not observed on 

the rocket flight; there is a possibility that it may eventually be obtained 

from satellite data. However, in order to continue the discussion we make 

the assumption that the spectrum is exponential, typical of soft auroral 

radiation, having the form. J(>E) = J o exp(-E/E ). Then, using the flux 

between 2 and 10 keV, deduced from the ionization rate, and the measured 

flux (>70 keV) of 1 x 10 3 cm 2 sec 1 ster" 1 , we find that J = 4 x 10^ cm -2 

-1 -1 

sec ster and E =8.3 keV, for E in units of keV. This value of E is 

° o 

within the limits of other observations of auroral spectra in the energy 

range, thus giving confidence to the assumed form of the spectrum. 

The specification of the spectrum allows the flux of electrons (>40 keV) 

at the time of our rocket flight, with K = 8, to be estimated as 
4 -2 -1 -1 

3 x 10 cm sec ster . This may be compared with a rocket observation by 
Gough and Collin (1973) at South Uist, geomagnetic latitude 58°: on one 
occasion, with = 8+, they found the flux (>40 keV) to be 1 x 10 4 cm" 2 sec" 1 
ster 1 . The difference between these two values is small considering the 
variability of the flux. Thus the agreement is considered satisfactory. 
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S.3 Ionization rate on other nights 

We believe , as was just discussed, that the ionization source for the 
upper E region on the very disturbed night can be identified as energetic 
electrons. It is interesting to consider next the extent to which energetic 
electrons contribute to the ionization on other, quieter, nights. There are 
now five observations of the intermediate layer at Wallops Island near mid- 
night: 12 April 1963 (Cartwright, 1964); 22 June 1965 and 22 February 1968 . 
(Smith, 1970); and 11 September 1969 and 1 November 1972, shown here in 
Figure 2. 

Applying the method of analysis described in section 4.3, we have 
derived the ionization rate for each of these five cases, for the lower and 
upper parts of the layer, and for the whole layer. These values are given 
in Table 1, which also includes the altitudes of the valleys and the peaks 
used as limits in the integration of the continuity equation, and the 

values of K for the three-hour period including the launch time. 

P 

Note the relatively small difference for a given night in values of 
ionization rate in the lower and upper part of the region. This implies that 
the ionization rate is nearly constant, at midnight, over the altitude range 
of 120 to 180 km. There is, however, from night to night, a variation of . 
almost three orders of magnitude. 

These ionization rates are plotted in Figure 4 against the corresponding 
value of K . The larger circles represent data from the whole layer; those of 
the lower and upper part are represented by smaller circles . In one case all 
three values coincide. The ionization rate in the intermediate layer at 
midnight is clearly related to the level of geomagnetic activity as represented 
by the least-squares fit: 

-3 -1 

q a 0.14 exp (0.71 K ) cm sec . (7) 

The correlation coefficient for the five observations is 0.96. 
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If the ionization is assumed to be caused by energetic electrons, equation 
(7) may be interpreted as descriptive of the variation with geomagnetic 
activity of the flux of electrons (>2 keV) at Wallops Island: 


J(>2 keV) M x 10 


exp (0.71 K ) 


-2 -1 -1 
cm sec ster 


( 8 ) 


A direct test of this relation may be possible from existing satellite 
data; it could certainly be tested in future rocket experiments. For the 
present, however, we can only suggest that it is reasonable, based on compari- 
son with a series of observations by Gough and Collin (1973). They measured 
the flux of electrons (>40 keV) at South Uist on 13 rocket flights and found 
a good correlation with geomagnetic activity. Their measured flux is related 
to at launch time by the least-squares fit: 


J (>40 keV) = 5.0 exp (0.98 K ) 


-2 -1 . -1 
cm sec ster 


( 9 ) 


With a correlation coefficient of 0.79. . 

Allowing for some variation of the energy spectrum from one occasion to 
another, the flux of electrons (>2 keV) and the flux of electrons (>40 keV) 
at South Uist show similar dependence on K . This supports our assumption 
that energetic electrons are the principal source of ionization in the upper 
E region near midnight at midlatitudes. 

5.4 Latitude variation 

Rowe (1973) has presented data from Arecibo, Puerto Rico, geomagnetic 
latitude 30°N, which show that the electron density of the upper E region 
at night is rather well correlated with geomagnetic activity. In this section 
we consider the relation between his measurements and the rocket observations 
at Wallops Island. 

The intermediate layer, seen regularly at Wallops Island, is not a 
feature of the nighttime ionosphere at Arecibo, presumably because of the 
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difference in wind systems at the two latitudes. The electron density at 
Are ci bo is nearly constant in the region 130 to 160 km. We have calculated 
the ionization rate there from electron density data using a value of recom- 
bination coefficient which is a function of altitude and K , as described in 
section 5.1. We shall consider first a simultaneous observation of electron 
density at Wallops Island and Arecibo. 

At our suggestion observations were made at Arecibo on the very disturbed 
night, 31 October - 1 November 1972. At the time of the rocket observation at 
Wallops Island, the electron density between 140 and 160 km at Arecibo was 
constant with a value of 6 x 10~ cm 2 (Rowe, personal communication). Hie 

average recombination coefficient over this range of altitude for K = 8 

P ' 

is calculated to be 1.8 x 1CT 7 cm 3 sec' 1 so that, for a steady-state condition, 
the ionization rate is 6.5 cm' 3 sec' 1 . This may be compared with the value 

_3 _i 

35 cm sec obtained from simultaneous observations at Wallops Island. The 
decrease of a factor of about 5 between Wallops Island and Arecibo may be 
attributed to the change in latitude. 

The ionization rates at Arecibo for 17 nights with various values of K 

P 

have been calculated from the average value of electron density between 130 and 
160 km, given by Rowe (1973), again using a value of recombination coefficient 

which includes the variation of electron temperature with altitude and K We 

P* 

find: 

q = 2.0 x 10“ 2 exp (0 . 77Kp) cm ^ sec (10) 

with a correlation coefficient of 0.79. 

Comparison of the variation of ionization rate with K at Arecibo, as 
given by equation (10) with the corresponding equation (7) for Wallops Island 
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reveals a remarkably parallel dependence on K . This is surprising in view of 
the relatively low geomagnetic latitude of Arecibo. However, it does suggest 
that electron precipitation may be an important source of ionization at night 
even on the equator-ward side of the sub-auroral zone. 

Observations at low latitudes of electrons with energies of a few keV 

have been reported by Hill et al . (1970) and by Hayakawa et al . (1973). The 

former, from rocket flights at Kauai, geomagnetic latitude 22°N, on two quiet 

days (K = 1+ and 0+) , find the flux of electrons with energies between 1 and 

-2 -1 -1 

10 keV to be about 100 cm sec ster at an altitude of 150 km. This 
flux appears to be smaller, by an order of magnitude, than the flux required 
to produce the ionization rate deduced from the Arecibo observations, for 
these very quiet conditions. 

6. Conclusion 

Analysis of electron density profiles in the upper E region near 
midnight at Wallops Island has indicated that the ionization rate is very strongly 
correlated with geomagnetic activity. We have interpreted this as evidence 
that the principle source of ionization under these conditions is energetic 
elections on geomagnetically quiet nights as well as on disturbed nights. 

Electron density data from Arecibo indicates a similar correlation between 
ionization rate and geomagnetic activity. Evidence of sufficient particle 
precipitation to account for the deduced ionization rates at Arecibo is 
lacking, however. The entire subject of precipitation of energetic electrons 
at latitudes below the sub-auroral zone is not resolved and requires further 
investigation. 

No additional source of ionization from scattered solar UV radiation is 
needed to account for the ionization rates at Wallops Island. In fact, the 
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present data indicate that an upper limit for ionization sources other than 

-1 -3 -1 

energetic electrons is 10 cm sec in the upper E region near midnight. 

The observations at Wallops Island also lead to the conclusion that 
energetic electrons are negligible as a source of ionization in the upper E 
region in the daytime. Even under the most disturbed conditions, energetic 
electrons contribute of the order of one percent of the daytime ionization rate. 
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FIGURE CAPTIONS 


Figure 1. Count rate observed on a disturbed night at Wallops Island using 
a Geiger counter sensitive to electrons (>70 keV) . 

Figure 2. Electron density profiles at Wallops Island near midnight on quiet 
and disturbed nights. 

Figure 3. The effect of winds in the E region. Calculated for a vertical 

ion velocity given by w = 20 sin[7r(150 - z)/30] m sec~^ for alti- 

-7 3 

tude z in km, with the recombination coefficient 1.9 x 10 cm 

sec \ for production rates of (1) 0.4, (2) 4.0, (3) 40, (4) 400 
-3-1 

and (5) 4000 cm sec . The thin line is the solution neglecting 

diffusion; the thick line is the complete solution with a diffusion 

8 2 - 1 

coefficient of 2 x 10 cm sec 


Figure 4. The variation with of the ionization rate in the upper E region 
from observations of the intermediate layer near midnight. 
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Table 1: 


Ionization Rate in the Upper E Region near Midnight. 


Altitude (km) 



Date 

Valley 

Peak 

Valley 

12 

Apr. 

1963 

120 

140.5 

182 

22 

Jun . 

1965 

125 

153 

177 

22 

Feb. 

1968 

116 

140 

154 

11 

Sep. 

1969 

118 

148.5 

180 

1 

Nov. 

1972 

115 

143 

175 


Ioni zation 

Rate (cm 

-3 -1, 

sec ) 


Lower 

Upper 

Whole 

K 

P 

0.90 

0.48 

0.62 

1 

0.090 

0.071 

0.081 

0+ 

2.0 

3.3 

2.5 

3+ 

0.60 

0.59 

0.59 

2 + 

26 

43 

35 
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